import numpy as np
import scipy.stats as stats
import scipy.optimize as opt